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The present study brings forward important information, within the framework of spectral distri- 
bution theory, about the types of forces that dominate three realistic interactions, CD-Bonn, CD- 
Bonn+3terms and GXPF1, in nuclei and their ability to account for many-particle effects such as the 
formation of correlated nucleon pairs and enhanced quadrupole collective modes. Like-particle and 
proton-neutron isovector pairing correlations are described microscopically by a model interaction 
with Sp(4) dynamical symmetry, which is extended to include an additional quadrupole-quadrupole 
interaction. The analysis of the results for the I/7/2 level shows that both CD-Bonn+3terms and 
GXPF1 exhibit a well-developed pairing character compared to CD-Bonn, while the latter appears 
to build up more (less) rotational isovector T = 1 (isoscalar T — 0) collective features. Furthermore, 
the three realistic interactions are in general found to correlate strongly with the pairing+quadrupole 
model interaction, especially for the highest possible isospin group of states where the model inter- 
action can be used to provide a reasonable description of the corresponding energy spectra. 



I. INTRODUCTION 

A unified microscopic description of light to heavy nu- 
clei requires a comprehensive understanding of the strong 
interaction and how it manifests itself in the nuclear 
medium. Effective interaction theory attempts to model 
the essence of this strong interaction in terms of one, 
two, and sometimes higher (three or even four)-body 
interactions for the purpose of supporting microscopic 
shell-model calculations that target reproducing striking 
features of nuclei. Such are strong pairing correlations 
found near closed shells that yield to collective rotational 
motion as one moves away from shell closure as well as 
more subtle effects that must be understood at a deeper- 
level to reproduce, for example, nuclear abundances as 
realized through fast or slow decay processes within and 
between nuclear species. While good progress is being 
made toward understanding the strong force, especially 
through the recent work on lattice quantum chromody- 
namics studies, much work remains to be done. 

Short such a comprehensive understanding of the 
strong interactions, one way to gain insight into the prin- 
cipal characteristics of various microscopic interactions is 
to perform a detailed study of their matrix elements or 
common quantities they generate, such as eigenvectors 
or eigenvalues. For instance, in the same environment 
two very similar interactions are expected to yield sim- 
ilarpatterns of their matrix elements (see, for example, 
[ESQ); as weu to yield large overlaps of corresponding 
eigenstates (e.g., 0, Q) and close energy spectra (e.g., 
0, 0)- On the other hand, a complementary compari- 
son that is based on the theory of spectral distributions 
[1, Q that invokes overall correlations of two interactions 
offers a broader view on their global behavior and uni- 
versal properties 0, El El, EMI • 

The theory of spectral distributions of French and col- 



laborators is an alternative approach for studying effec- 
tive interactions [1, El, EH and continues to be a powerful 
concept with recent applications in quantum chaos and 
nuclear astrophysics including studies on nuclear level 
densities, transition strength densities, and parity/time- 
reversal violation (for example, see [T3, El, El, l20l . [HI, 
HH). The significance of the method is related to the 
fact that low-order energy moments over a certain do- 
main of single-particle states, such as the energy centroid 
of an interaction (its average expectation value) and the 
deviation from that average, yield valuable information 
about the interaction that is of fundamental importance 
[HI, El, HI, HI, [H, |H, [23, [H, ^ without the need for car- 
rying out large-dimensional matrix diagonalization and 
with little to no limitations due to the dimensionality of 
the vector space. Within this general framework, a sim- 
ple and elegant global comparison of pairs of interactions 
can be carried in a unified way regardless of how the in- 
teractions are built or of the models that adopt them, 
ft is based on their correlation, which is a measure that 
is independent of the centroids of the interactions. This 
correlation measure turns on a comparison of the one- 
and two-body parts of the interactions, and in so doing 
probes beyond the mean-field potential. 

In addition, the group-theoretical foundation of the 
spectral distribution methods naturally establishes a 
propagation of information from nuclear systems of two 
particles to many-fermion nuclei [1, El] and makes the 
approach especially suitable for studies of the good- 
ness/breaking of symmetries in complex many-particle 
systems [H, El HH ■ Such studies can likewise help reduce 
the dimensionality of a model space to a tractable size 
by detecting the optimal subspace domain for a partic- 
ular many-body problem where microscopic calculations 
become feasible d, [H, HI] • 

In this paper we employ the theory of spectral distri- 
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butions to compare three realistic interactions, namely 
CD-Bonn [Hj], CD-Bonn+3terms [H and GXPF1_[| 
and two pairing+quadrupole model interactions [H, 1341 ]. 
Such a study is important for understanding the types 
of forces that dominate a realistic interaction and its 
ability to describe correlated and collective phenomena. 
Most significant are the formation of nucleon pairs and 
quadrupole collective excitations that possess a clear 
symplectic algebraic structure, sp(4) and sp(3,R), re- 
spectively. 

The Sp(4) dynamical symmetry of like-particle and 
proton-neutron pairing correlations (33l . l35l . l36l | between 
nucleons occupying the same major shell has been found 
to provide for a reasonable microscopic description of the 
pairing-governed isobaric analog + states in light and 
medium mass nuclei [341 . |37| . Currently, these nuclei have 
a significant impact in astrophysical studies [H, [39[ and 
have been extensively explored after the advent of ra- 
dioactive beams. The comparison with realistic inter- 
actions can determine the extend to which the signifi- 
cantly simpler Sp(4) model Hamiltonian can readily be 
used to obtain an approximate, yet very good description 
of low-lying nuclear structure and in turn, one can ap- 
ply the model to larger model spaces that are otherwise 
prohibitive in size. Furthermore, we introduce a possible 
Sp(4) symmetry breaking by an additional quadrupole- 
quadrupole interaction and examine the capacity of the 
extended model interaction to imitate realistic interac- 
tions. This, in turn, provides a further step towards gain- 
ing a better understanding of the underlying foundation 
of the microscopic interactions. 



II. SYMPLECTIC sp(4) PAIRING MODEL 
INTERACTION 

The close interplay of like-particle and proton-neutron 
isovector (isospin =1) pairing correlations have long been 
recognized as a major driver that shapes nuclear systems 
with valence protons and neutrons occupying the same 
major shell. While like-particle pairing interactions are 
known to dominate far from the N = Z line, closer to 
it strong proton-neutron pair correlations are also very 
important. Hence, isovector pairing plays a crucial role 
in understanding the microscopic structure of light and 
medium mass nuclei around as well as far off the val- 
ley of stability. A group-theoretical microscopic descrip- 
tion of isovector pairing, based on the fermion realiza- 
tion of the so (5) algebra [35| (isomorphic to sp(4)), was 
successfully applied to the structure of /p-shell N = Z 
nuclei [36]. These algebraic results have since been con- 
firmed through pairing-plus-quadrupole shell model work 
[4(i| . Indeed, some recent results show that the symplec- 
tic Sp(4) dynamical symmetry is fundamental to the nu- 
clear interaction that governs fully-paired isobaric analog 
+ states of light and medium mass even- A nuclei with 
valence protons and neutrons occupying the same shell 
0. 



The general model Hamiltonian with Sp(4) dynamical 
symmetry for a system of n valence nucleons in a 4f2- 
dimensional space consists of one- and two-body terms 
and can be expressed through the Sp(4) group generators, 

H sm = -GEL-i A\A t - FAlAo - £(f* - M) 

-D(fS-^)-C^l-eN, (1) 

where N counts the total number of valence particles, 
f' 2 = Q,{f+,f^} + T 2 is the isospin operator, ^J.+i.-i 
creates a proton-neutron (pn) pair, a proton-proton (jyp) 
pair or a neutron-neutron (nn) pair of total angular mo- 
mentum J w = + and isospin T = 1, G,F,E,D and 
C are interaction strength parameters and e > is the 
Fermi level energy. This Hamiltonian, which is rotation- 
ally invariant, conserves the number of particles and the 
third projection (To) of the isospin, while it includes scat- 
tering of a pp pair and a nn pair into two pn pairs and 
vice versa, along with a J-independent isoscalar (T = 0) 
pn force. The significant interplay between isovector and 
isoscalar interactions is evident in the low-lying struc- 
ture of N = Z odd-odd nuclei with valence protons and 
neutrons filling the same major shell. 

Estimates for the interaction strength parameters in 
(fl"j) were found (33l . H3| as a result from an optimal re- 
production of the Coulomb corrected [4l| experimental 
energies [U |43| of the lowest isobaric analog + states of 
even- A nuclei with valence nucleons occupying the I/7/2 
orbit or the l/ 5 /22pi/ 2 2p3/2lg9/2 major shell [55|. For 
the I/7/2 level with a 40 Ca core the interaction strengths 
were estimated to be, G/Q = 0.453, F/Q = 0.072, C = 
0.473, D = 0.149, E/{2Q) = -1.120, e = 9.359. The 
analysis revealed that the model interaction with Sp(4) 
dynamical symmetry accounts quite well for the available 
experimental energies of isobaric analog + states for a 
total of 149 nuclei [34[ and in addition for the observed 
detailed structure beyond mean-field effects such as the 
N = Z anomalies, isovector pairing gaps and staggering 
effects [37| • This in turn allowed us to interpret the main 
driving force that defines the properties of the states un- 
der consideration and to provide a reasonable description 
of these states, while retaining the physical validity and 
the proper limits of the strengths of interactions available 
in literature. 

An important feature of our algebraic Hamiltonian ([I]) 
is that it arises naturally within a microscopic picture. 
Because of this, the Sp(4) interaction can be compared 
to realistic interactions and, as well, the physical nature 
of the model interaction and its strength can be real- 
ized. From a microscopic perspective, the pair-creation 
operators, A^ , and their annihilation counter parts, A, 
are realized in terms of creation Cj mrT and annihilation 
Cjma single-fermion operators with the standard anti- 
commutation relations {c jmcr , cj, mV ,} = tfjj'tfm.ro'^a-.o-'' 
where these operators create (annihilate) a particle of 
type a = ±1/2 (proton/neutron) in a state of total an- 
gular momentum j (half integer) with projection m in a 
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finite space 2f2 = Sj (2j + 1) . There are ten independent 



scalar products (zero total angular momentum) 
fermion operators: 



of the 
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A„ 



T+ 



N 



1 



^20(1 + ^0 



(4) f . G« 



i 



J"l" J, 



/2fi 



^ C ]™,±l/2 C i m 'Tl/2' 



= E*E 



(2) 



which form a fermion realization of the symplcctic Bp (4) 
Lie algebra. Such an algebraic structure is exactly the 
one needed to describe isovector (like-particle plus pri) 
pairing correlations and isospin symmetry in nuclear iso- 
baric analog + states. 

Using relations j2]), the one- and two-body interaction 
([!]) can be rewritten in standard second quantized form in 



terms of fermion creation 



nihilation a jm ^/2)a = (-1) 
operators, 



l jm(l/2)a 
i+m+l/2+u 



>(l/2)" 



and an- 



Cj- m (i/2)-o tensor 



H = — 



t < u 

r = (j, T) 

E VF^rM ® « s } ( ° 0) - 7 E V(l + ^)(l + <5 tu )[r]W r r atu {{at ® at} r ® { at ® au } r }( 00 ), (3) 



r < s 
{>~=(jV,|)} 



rstu 

r 



with [r] = 2(2j' r + l) and [r] = (2J+1)(2T+1), where e r 
is the single-particle energy and W/£ u * s the two-body 
antisymmetric matrix clement in the JT-couplcd scheme 
.(_)r+.-r w r 



srtu 

{—) r+s ~ t ~ u Wg rut = W[ urs \. For an isospin nonconserv- 
ing two-body interaction of isospin rank T, the coupling 
of fermion operators is as follows, {{a£ <E> o) s } JT ® {a t <8> 
a«} JT } (0T) , with V^^ 77, matrix elements. The latter 
are expressed through the parameters of the model in- 
teraction for isospin rank and 2 of H S pu\ and {r < 
(s,t); t < u} orbits as follows, 



w. 



(0)JT _ 



WJL = < rsJTMT () \H (0) \tuJTMTn > 



W. 



(2)JT 



-(G + -gO ^—^$(JT),(01)3rsdtu 

H ^ + i )[HT + ^ 1 + c}u " (4) 

< rsJTMT \HW\tuJTMT > 

= ^-(F^—^—^-Sj S rs Stu — DS r tS su )S T i- (5) 

The isotensor part ([5|) of the model interaction intro- 
duces isospin dynamical symmetry through the _D-tcrm 
(retaining T as a good quantum number and splitting 
the energy degeneracy along the third projection of the 
isospin) and as well a plausible, but very weak, isospin 
mixing (F-term) |44j . 

For the purposes of this paper, we will use only the 
isoscalar part of our model Hamiltonian (0| and set all 
the orbits equal to j = 7/2 (r = s = t = u) because we 
choose to focus on a study of nuclei in the single 1/7/2 



level. In addition, these matrix elements correspond to 
the pure nuclear interaction and do not include Coulomb 
repulsion because its effect is corrected in the experimen- 
tal energies themselves by applying an empirical formula 
deduced in 41|. This may result in slightly more bound 
states predicted by our model when compared to esti- 
mates of realistic interactions. 

Within the isospin-invariant picture, the two-body ma- 
trix elements of the model Hamiltonian W^ u ^ depend 
only on three parameters, Go = G + Eq = + 
and G, 



7 7 



-Gq8(jt),(o 



^-{-EoK-f + ^+C}. (6) 



The two-body matrix elements reflect the microscopic as- 
pect of the model interaction, which is J-indcpendent for 
all but J = 0. Hence H sp ^ describes the average be- 
havior of higher- J states, while it distinguishes between 
T = and T = 1 groups of states. The smaller the mag- 
nitude of Eq (< 0), the smaller the separation of these 
groups. As expected, the pairing correlations contribute 
only to the (J = 0, T — 1) state and they are absent 
for higher-J states where both particles are uncoupled. 
Relative to the + T = 1 state, the bigger the Go pair- 
ing strength, the larger the energy gap to the highcr-J 
states. 

The role of the Sp(4) dynamical symmetry in generat- 
ing the energy spectrum of the I/7/2 nuclei can be further 
understood by comparing the S p(4 ) interaction to the 
CD-Bonn[Hl, CD-Bonn+3terms M and GXPF1 re- 
alistic interactions. CD-Bonn is a charge-dependent one- 
boson-exchange nucleon-nucleon (NN) potential that is 
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one of the most accurate in reproducing the world proton- 
proton and neutron-proton scattering data. In addi- 
tion, the CD-Bonn+3terms interaction introduces phe- 
nomenological isospin-dependent central terms plus a 
tensor force with strengths and ranges determined in no- 
core 0Hlu shell model calculations to achieve an improved 
description of the A — 48 Ca, Sc and Ti isobars. The 
GXPF1 effective interaction is obtained from a realis- 
tic G-matrix interaction based on the Bonn-C potential 
[5(j by adding empirical corrections determined through 
systematic fitting to experimental energy data in the fp 
shell. 



III. THEORY OF SPECTRAL DISTRIBUTIONS 

Group theory underpins spectral distribution theory 
d, [H, [45[. The model space is partitioned ac- 

cording to particular group symmetries and each sub- 
sequent subgroup partitioning yields finer and more de- 
tailed spectral estimates. For n particles distributed over 
Af single-particle states, a scalar distribution (denoted 
by "n" in the formulae) is called the spectral distribu- 
tion averaged over all n-particle states associated with 
the U(7V = 4£7) group structure and an isospin-scalar 
distribution (denoted by "n,T") is averaged over the en- 



semble of all n-particle states of isospin T associated with 
U(A/" = 20)<X>U(2) T . 

For a spectral distribution a (a is n or n,T), the cor- 
relation coefficient between two Hamiltonian operators, 
H and H' , is defined as 

<=H,H' ~ 1 — Z \'J 



(H^H') a - {H^) a {H') a 
(Jhvh> 



(8) 



where the "width" of the distribution is the positive 
square root of the variance, 

(a£) 2 = ((H (H) a )Y = (H 2 ) a ((H) a ) 2 . (9) 

The average values, related to the trace of an operator 
divided by the dimensionality of the space, are given in 
terms of the ensemble considered. In the (isospin-)scalar 
case, the correlation will be denoted by £™ (C"' T ) or sim- 
ply C (C T ) for n = 2. 

The steps for computing the C" correlation coefficient 
and the a a variance (56j are given in 
(see also computational codes [ID, |46[) and take on the 
simple form for a single-j level: 



(H^HT-(HY(HT =E^(«,r)^^^[A]IT/ r ; r (2)IT' r J ; rr (2), (10) 

w/ r ; r (2) =w/ r T rr -wy\ (ii) 

=^ fXT ^[A]W^ rr , (12) 

where r = {0 or 1} is the isospin label of the two-body matrix elements, W^ 7 T rr (2) is the traceless pure two-body 

interaction and is the monopole moment or centroid in the (isospin-)scalar case, that is the average expectation 
value of (the isospin-r part of) the two-body interaction for a two-particle system n = 2. In the scalar (a = n) case 
the following holds, A = T = (J,t), = CsD' ^ = 40 = 2(2j r + 1), and the r- independent propagator is 

n(n-l){N-n){N-n-l) 
^ T) = 2(Af-2)(Af-3) • (13) 

In the isospin-scalar (a — n, T) case: A — J, ^2j[J] = K^+izllA^ j\f = 2S1, and the propagator functions are [13, S3 
M n,T,r = 0) = ["(" + 2)-4T(r + l)][^^-^l)-nT + l)] (M) 

+jjn(n -2) { N-lW- n - + 1){N + 1){N + 2) 

+±T(T + 1)[(57V — 3)CAT + 2)n(f - N) + M(M - + + 6)]}. (15) 

I 



In terms of a geometrical picture, the correlation coeffi- 



cient ( defines the angle between two vectors (H and H') 
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of length a H (,) ([9]) and hence its square gives a normalized 
measure (percentage) of one of the vectors, e.g. the Sp(4) 
interaction, that is contained in the other, such as a real- 
istic interaction. The correlation coefficient is a measure 
that is independent of the averages of the interactions. 
Clearly these averages, though an interesting measure, 
are irrelevant when the focus is on detailed property- 
defining two-body interaction beyond strong mean-field 
effects. 

For the Sp(4) interaction, the average two-body inter- 
action is expressed in terms of the model parameters in 
the scalar case as, 



Wc 



3En 



3 Go 

+ 2(^-1) 



C 



(16) 



and in the isospin-scalar case as 
Go 



S T1 +E [(-iy +±]-C. (17) 



Hence the pure two-body Wtj 7 7 7 (2) matrix elements 

2 2 2 2 

(|6"|), and consequently the correlation coefficients involv- 
ing i? s p(4), are independent of the C (and Eq) parame- 
ter(s) in the (isospin-)scalar case. 



IV. UNDERLYING SYMMETRIES OF 
REALISTIC INTERACTIONS 

We now use statistical concepts to pro be the nature of 
the CD-Bonn [p], CD-Bonn+3terms [H and GXPF1 @ 
realistic interactions, hereafter referred as Hr. Specifi- 
cally, we will compare these interactions to the symplcc- 
tic pairing and quadrupole interactions through their mu- 
tual correlations. Clearly, if two interactions have similar 
matrix elements they will be strongly correlated and any 



pattern that is observed in the behavior of one will be 
reflected in the other. This can be made quantitative by 
evoking measures from statistical spectroscopy, namely, 
the closer the correlation coefficient between two interac- 
tions is to unity the more similar their spectra with the 
two coinciding within a rescaling factor when the corre- 
lation coefficient is unity. 

In a similar manner, the projection onto a model 
Hamiltonian that describes collective rotational excita- 
tions or/and pairing correlations can be used to probe 
the rotational and pairing characteristics of a microscopic 
interaction [ll|, HE [H, |49[ . The dynamical symmetry of 
the pairing (or quadrupole-quadrupole) interaction sets 
a specific relation between the matrix elements of the 
Hamiltonian that models it. If this relation is found in a 
realistic interaction, that is, the model and realistic inter- 
actions are strongly correlated, then the latter possesses 
the underlying symmetry and will reflect the characteris- 
tic properties of the pairing (quadrupole) Hamiltonian. It 
should be clear that the complement is also true, namely, 
if a model interaction is strongly correlated with a real- 
istic one, the associated model calculations can be used 
to investigate the behavior of physical systems. 



A. The Sp(4) Model and Pairing Character 

An interesting feature of any interaction is its trace- 
equivalent part. If the latter is found dominant then only 
the underlying group scalars are enough to provide for an 
approximate and yet reasonably good solution The 
greatest advantage in this case is the simplicity of the 
many-body problem and the tractable size of the model 
space. In the isospin-scalar case, the centroid of a Hamil- 
tonian expressed through the e single-particle energy and 
the W®' 1 monopole moments (fT2"|) is [47[ 



J 



-ne 



- [T(T+ 1) 



r 



3 ,W}-W° 



n 

4 J 



(18) 



For a Hamiltonian with symplcctic dynamical symmetry, 
H S p(4), the trace-equivalent part in the isospin-scalar dis- 
tribution includes the E-, C- and e-terms of ([1]). When 
applied to the lowest isobaric analog + states of the nu- 
clei in the I/7/2 orbit, it reproduces their energy within 
1% of the experimental value for about a third of the 
nuclei. While for these states the centroid is sufficient 
to achieve a good description, its difference with experi- 
ment goes up to 7% compared to only 0.4% achieved by 
the whole Sp(4) Hamiltonian. In addition, a model with 
a trace-equivalent Sp(4) Hamiltonian will not be capable 
of explaining the fine nuclear structure where H B pu\ (p} 



succeeded |37( and will not correlate with any of the real- 
istic interactions. The latter indicates an inadequate re- 
production of the entire energy spectrum. Indeed, while 
such an interaction was found insufficient for a descrip- 
tion of ds shell nuclei when compared to several effective 
interactions, a drastic improvement was achieved with 
the inclusion of pure two-body residual interactions of 
the spin-orbit and quadrupole types [ll[ as well as pair- 
ing correlations (26[. In summary, the Sp(4) symmetric 
Hamiltonian |T]) provides for a more accurate description 
of nuclear structure by adding to an average interaction 
suitable for the isobaric analog + states in I/7/2 a sig- 
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nificant isovector pairing part. 

Furthermore, the Sp(4) dynamical symmetry allows 
the model Hamiltonian to reflect on the charge de- 
pendence of the nuclear interaction, which is evident 
from experiments and present in almost all of the mod- 
ern realistic interactions (e.g., CD-Bonn). While the 
small isospin admixture found in the I/7/2 isobaric ana- 
log + states has been directly estimated through the 

sp(4) eigenstates 44 1, the theory of spectral distribu- 



tions provides a further estimate of isospin symmetry 
breaking throughout the entire spectrum (4?| based on 
the Sp(4) isotensor interaction ([5]). Using equations (4), 
(15) and (25) in [l(j the isospin T + 2 admixture into 
an average T state is found in I/7/2 to be, as expected, 
much smaller (on average less than 0.0001%) than the 
one detected among the + seniority-zero states [H| . As 
expected, it is also much smaller than the measure calcu- 
lated in for the two-body Coulomb interaction in the 
I/7/2 shell (with a maximum value of 0.009%), because 
the latter corresponds to the stronger AT = 1 admix- 
ture. The quite small isospin symmetry breaking that 
the isotensor Sp(4) model interaction introduces allows 
us to carry the present study without its consideration. 

The extent to which the Sp(4) dynamical symmetry 
governs the Hr realistic interactions within a certain do- 
main of states is represented by the correlation coeffi- 
cients between H 5p ^ and Hr (Table UJ. 



TABLE I: Correlation coefficients for a two-nucleon system, 
n — 2, in the scalar (£) and isospin-scalar (£ T ) distributions. 

C CD-Bonn CD-Bonn+3terms GXPF1 

(C T=0 ,C T=1 ) 

#sp(4) 

CD-Bonn 
CD-Bonn+3terms 



0.66 
( - ,0.61) 



0.64 
( - ,0.85) 

0.95 
(0.99,0.94) 



#q(2) 
Hm 



0.47 
(0.60 ,0.73) 

0.81 
(0.60 ,0.95) 



0.60 
(0.68 ,0.50) 

0.87 
(0.68 ,0.98) 



0.76 
( - ,0.71) 

0.96 
(0.98,0.99) 

0.97 
(0.99,0.97) 

0.53 
(0.74,0.65) 

0.93 
(0.74,0.96) 



In the scalar distribution, where the outcome is aver- 
aged over the isospin values, the analysis of the results 
shows that all of the realistic interactions correlate be- 
tween themselves to a high degree (Table |H upper cell 
for each pair of interactions). In comparison, each of 
them has a correlation with the Sp(4) symmetric inter- 
action of order of 0.6-0.8, which is typically regarded as 
a good one [HI]. This implies that the realistic interac- 
tions possess around 40-60% of the dynamical symplectic 
symmetry of # 5p(4 ) [((H sp{ih H R ) 2 } and hence 0.4-0.6 por- 
tion of Hr is dynamically symmetric under Sp(4) trans- 



sideration. This is a very interesting result, and definitely 
valuable concerning the restrictions the symplectic model 
is subject to. 

A much more interesting scenario occurs when the 
isospin-scalar case is considered. This is because in 
this case the space is divided into two regions specified 
by their isospin values with a view towards a more de- 
tailed examination of the nature of the interactions under 
consideration. Indeed the centroids of both T = and 
T = 1 regions are considerably separate as is observed in 
the realistic and Sp(4) interactions and as well confirmed 
by experiment. In addition, the important pairing corre- 
lations that are described in the symplectic model enter 
in the T = 1 channel, where the tendency towards pair 
formation of realistic interactions can be detected. 

As in the scalar case, all of the realistic interactions are 
quite strongly correlated in both the T = and T = 1 
channels (Tablcfl] lower cell for each pair of models). The 
T = 1 correlation coefficients between £f S p(4) and the Hr 
realistic interactions do not depend on any of the param- 
eters in H B pu\ including the pairing strength itself [57j | 
and show that H B p(4\ correlates strongly with the realis- 
tic interactions. Among the three Hr, the T = 1 part of 
the CD-Bonn+3tcrms interaction possesses the closest 
similarity to the pn and like-particle J — pairing corre- 
lations. This is indicated by its large projection of 72% 



in Table H] onto the T = 1 H sp / 4 \ pairing 



formations. Equivalcntly, the if. 



0p(4) 



model interaction 



contains 40-60% of the realistic interactions under con- 



interaction. Hence, the CD-Bonn+3terms interaction is 
expected to describe quite well phenomena of a pairing 
character. 

The individual pairing strength associated with each 
realistic interaction is typically invoked for purposes of 
comparison. Compared to CD-Bonn, the J = isovec- 
tor pairing strength estimate turns out to be stronger 
for both CD-Bonn+3terms and GXPF1 with a rela- 
tively weaker coupling observed in the latter. Indeed, 
GXPF1 was shown to tend towards smaller J = pairing 
strength [2j when compared to two other effective inter- 
actions, namely the G interaction based on the Bonn-C 
potential [5(| and the KB3G interaction [5l|. However, 
pairing effects, with strong or weak coupling, may be 
fully or partially reflected in Hr considerably depending 
on the strength of the overall interaction. It is the corre- 
lation coefficient between a realistic (Hr) and symmetry- 
holding (as H B pu\) interactions that manifests what part 
of Hr is ruled by the symmetry. In this sense, we can 
identify that pairing features arc more fully developed in 
CD-Bonn+3terms, then in GXPF1 and the least in CD- 
Bonn ((jj Hr in Table |T|. In the latter, other types 
of interaction compete stronger with pair formation than 
in the former two interactions and hence suppress pairing 
coherence. 

The large J = coherence and its strong coupling ob- 
served in the CD-Bonn+3terms interaction should not 
be surprising because it reproduces (by an optimal fit) 
the energy difference between the ground state and the 
first 2 + state of 48 Ca. Such an observable is believed to 
be directly affected by the formation of correlated pairs 
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in the ground state of the spherical core of 8 Ca and the 
pairing gap that occurs below the first excited state of 
a broken pair. It is interesting to point out that while 
the addition of three phenomenological terms to the CD- 
Bonn interaction to obtain CD-Bonn+3terms keeps the 
close similarity between both interactions [C 1 = 0.94 (Ta- 
ble [I])] . it causes the correlations with isovector pairing 
interaction to double in strength 4 # H ) 2 ]- In short, 

the analysis shows that the simple Sp(4) model interac- 
tion can reproduce reasonably the T = 1 low-lying energy 
spectra generated by the CD-Bonn+3terms realistic in- 
teraction for a system of two nuclcons in the I/7/2 orbit 
and for this reason can be used as a good approximation. 

Another result in favor of the algebraic sp(4) model 
follows from a comparison of the lowest isobaric analog 
+ states in the A = 42 isobars. These are precisely the 
(J = 0, T = 1) states, which are expected to be shaped 
by stro ng p roton-neutron and like-particle pairing corre- 
lations |52| and are well described by the sp(4) model [HI]. 
The outcome reveals a very close similarity between the 
estimate of the two-body (J = 0, T = 1) matrix element 
for the symplectic interaction (—1.85 MeV) and both CD- 
Bonn+3terms (-2.06 MeV) and GXPF1 (-2.44 MeV) 
realistic interactions. In addition, the energy differences 
between the first 2 + state and the + ground state for the 
different effective interactions are also very close, namely, 
1.91 MeV (for H sp{i) ), 2.00 MeV (for CD-Bonn+3tcrms) 
and 1.50 MeV (for GXPF1). All these estimates are 
rather different from CD-Bonn with 0.48 MeV 2+ to 
+ energy difference but very close to the experimental 
energy gap for the A = 42 isobars, namely, 1.56 MeV for 
42 Ti, 1.59 MeV for 42 Sc and 1.52 MeV for 42 Ca. 



B. Pairing + Quadrupole Model Interaction 

While the pairing-governed isobaric analog + 
state energies are well determined within the framework 
of the Sp(4) model, the nuclear spectrum as described 
by the Sp(4)-symmctric Hamiltonian contains degener- 
ate higher- J states averaged for a given isospin value as 
can be clearly seen from its microscopic structure ((4]). 
Nonetheless, the correlation of -ff 5p (4) with realistic in- 
teractions for the I/7/2 level turns out to be reasonably 
strong. A question one can pose concerns the role of other 
significant interactions in nuclei such as the quadrupolc- 
quadrupole interaction (Q ■ Q). As our results indicate, 
answers to such questions can be found within the frame- 
work of statistical measures. 

The pairing model based on the sp(4) algebra [incor- 
porating like-particles pairing through an su(2) subalgc- 
bra] is commonly considered to be inappropriate for two 
reasons. The first reason is related to the degeneracy 
of the single-particle levels, which is not a problem for 
the I/7/2 shell considered as a single orbit well-separated 
from the ds shell and the upper fp shell. The second rea- 
son is the lack of the Q ■ Q interaction. This is because 
one usually neglects the fact that the pairing interaction 



contains in itself a part of the quadrupolc-quadrupole in- 
teraction. This part is not negligible with the correlation 
being typically between 0.4 to 0.6 depending on the dis- 
tribution considered: 15% when the whole space is con- 
sidered and 35% in the T = 1 region. This is probably 
one of the reason why the Sp(4) model interaction turns 
out to work rather well despite an explicit appearance of 
the quadrupolc-quadrupole interaction. 

Because of the fact that the Q ■ Q interaction is already 
present in the Bp (4) Hamiltonian, its additional influence 
can be studied following the construction prescribed in 
[ill ]. In short, we add a Q ■ Q term to the symplectic 
i?sp(4) Hamiltonian {1} in a way that this term is only 
the part of the pure two-body quadrupolc-quadrupole in- 
teraction that is not contained in Sp(4), or in the vector 
algebra terminology we add only the part that is orthog- 
onal to the pure two-body Sp(4) Hamiltonian [lq . 

H M = H 5P{4) + H^(2), H Q = -|Q • Q. (19) 

Such a Hamiltonian does not affect the centroid of 
Hsp(4) because Hq(2) is traceless. In this way this col- 
lective interaction preserves the shell structure that is 
built into -ff sp (4) and established by a harmonic oscilla- 
tor potential and as a result is favored in many studies 

0, tH, M, Hi. 

Compared to the pairing if 5 p(4) Hamiltonian, the ad- 
ditional collective interaction, Hq(2), has a lower cor- 
relation with Hr for all the cases except for T = 1 CD- 
Bonn and where Ch sp(4) ,Hr cannot be determined (Table 
U). The realistic interactions contain the Sp(4) interac- 
tion by 5% to 50% more than they contain Hq(2), the 
largest value being for the T = 1 CD-Bonn+3tcrms in- 
teraction. We should emphasize that this outcome does 
not imply that the realistic interactions correlate bet- 
ter with the pairing interaction than they do with the 
quadrupolc-quadrupole interaction nor that their pairing 
character is dominant. This is because the Hq{2) inter- 
action represents only that part of the rotational inter- 
action that is not included in the Sp(4) interaction and 

the entire Q • Q collective mode affects both clJ? „ 

(T) 

and Q H Hr correlations. The outcome only implies 
that a comparatively larger part of the overall correla- 
tions is already accounted for solely by the symplectic 
sp(4) algebraic model interaction. 

In our study, we vary only the quadrupole strength 
parameter in HH), to find its optimal value (which is an 
exact solution) by maximizing the correlation coefficient 
£ between Hm and Hr [54[. We do not alter the pa- 
rameters of the Sp(4) model, which have already been 
shown in an appropriate domain of states to be valid for 
reproducing various quantities (such as binding energies 
and pairing gaps) and are in agreement with estimates 
available in literature [U, H3] ■ This implies that the a 
"width" of -ff s p(4) © does not change. The minimization 
procedure is performed for Hm compared to each real- 
istic interaction and in the isospin-scalar case, for each 
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isospin value (Tables [T] and [II]) . 

TABLE II: First and second energy moments [the centroid W c 
(|12|) and the "width" a of realistic and model interactions 
for a two-nucleon system, n = 2, in the I/7/2 level. Hm 
is determined by an estimate for the quadrupole-quadrupole 
strength \ f° r each realistic interaction; its centroid energies 
coincide with the ones of H sf (4\ for a given distribution. 

CD-Bonn CD-Bonn+3tcrms GXPF1 

Scalar Distribution 



W c 


0.30 


0.19 


-0.60 


a 


0.59 


0.80 


0.98 


H M , Wo 


-0.63 


-0.63 


-0.63 


Hm, 


1.23 


1.36 


1.21 


Hm, X 


0.096 


0.124 


0.092 


Isospin-scalar Distribution, T — 


1 


wl 


0.54 


0.46 


-0.17 


a 1 


0.27 


0.57 


0.62 


H M , Wl 


-0.01 


-0.01 


-0.01 


Hm, a 1 


0.55 


0.41 


0.48 


Hm, X 


0.071 


0.036 


0.055 




(a) 



FIG. 1: (Color online) Geometrical representation of 
the realistic interactions, CD-Bonn (light blue), CD- 
Bonn+3terms (red) and GXPF1 (green), in an abstract op- 
erator space, where the horizontal plane is spanned by the 
orthogonal linear operators, the pure two-body H sv u\ and 
Hq(2) model Hamiltonians, both linearly independent of ad- 
ditional operators represented by the vertical axis, (a) Scalar 
distribution, (b) Isospin-scalar distribution, T = 1. The ori- 
entation of the vectors remains the same for any particle num- 
ber n > 2 in (a) and for all T = n/2 cases in (b). 



In both scalar and isospin-scalar cases, the addition 
of the quadrupole-quadrupole interaction definitely im- 
proves the (jy Hr correlation (Table [J) , which is asso- 
ciated with the angle between the two-body effective in- 
teraction and its projection on the plane spanned by two 
orthogonal vectors, two-body -ff S p(4) and Hq(2), in an 
abstract operator space (Figure [l}. In the same repre- 
sentation, the angles between the realistic interactions 
and both axes give the C^} TT and Cu u correla- 

(T) 

tions, and the length of each vector is specified by a H 
©. Therefore, while enhanced quadrupole effects rotate 
the projection of an interaction closer to the Hq(2) axis, 
greater influence of additional interactions neglected in 
Hm pushes Hr away from the Hm horizontal plane. 

The GXPF1 interaction is found to correlate best with 
the model Hm Hamiltonian for the scalar distribution 
(86%) and in the T = case compared to the other 
interactions (Table HJ). More than 50% of the T = 
GXPF1 interaction is accounted by the isoscalar model 
interaction. The T = correlation coefficients between 
the model Hm Hamiltonian and Hr are reasonably good 
and determined solely by the Q ■ Q interaction indepen- 
dent of its strength, x 033) • 

In the T = 1 region, all the realistic interactions con- 
sidered are reproduced to the 90% - 97% level by such 
a pairing+quadrupole model interaction. Other interac- 
tions contribute almost negligibly as is clearly seen from 
Figure [T] (b) with their contribution being least for the 
CD-Bonn+3terms interaction. In addition, within the 
T = 1 distribution, the smallest \ value (Table [IT} is 
found for the CD-Bonn+3terms realistic interaction as 
expected due to its reasonable correlation with the Sp(4) 
model interaction. In summary, the results once again 



prove that the pairing and quadrupole-quadrupole inter- 
actions are significant in shaping nuclear structure and 
are dominant for the T = 1 two-body nuclear interaction 
for the I/7/2 orbit. 

From the point of view of the model interaction (fT9^) 
we adopted, the rotational character of the three realistic 
interactions may appear to be obscure because a strong 
correlation to the entire Hq quadrupole-quadrupole in- 
teraction (not only to its projection Hq(2)), Ciiq h r i 
is needed. However, development of rotational features 

(T) 

turns out to follow qualitatively the C H ±^ 2 ) h corre l a_ 

tion coefficients (Tabic HJ). This is because Hq is already 
present in the Hm model Hamiltonian and one can vi- 
sualize Hq in the scalar (isospin-scalar, T = 1) case 
in Figure [1] as an axis that lies in the horizontal plane 
around 67° (54°) counterclockwise from the ff sp (4) axis 
and hence it lies furthest away from i? 6p (4) compared 
to all three realistic interactions. Therefore, the com- 
paratively largest collectivity is attained within CD- 
Bonn+3terms in the scalar case, within GXPF1 for 
the T — domain of states (with an exact quanti- 
tative measure C^ Q Hr = Cga^) h r )' an< ^ within CD- 
Bonn in the T = 1 case. In contrast, the T = 1 
part of the CD-Bonn interaction shows the smallest 
individual quadrupole strength, which is an illustrative 
example of a very prominent rotational behavior (de- 
tected via £) but of a weak strength (depending on a). 
In general, the individual quadrupole strength associ- 
ated with the GXPF1 interaction is very similar to, yet 
slightly stronger than, the one of CD-Bonn+3terms with 
the same trend observed for GXPF1 0] with respect to 
KBG3. 
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C. Energy spectrum tcraction with the realistic interactions and hence a sim- 

ilar pattern of energy states is expected. 

The scalar and isospin-scalar T = 1 distributions show 
strong correlations of the pairing+quadrupole model in- 



CD-Bonn 



> 



c 



5+.0 
~2+, 1 



-6+;? 



CD-Bonn+ 
— ] 3 terms 



-6+ 1 
-3+ 0j+ i 
■ 1+ ' 
= ' 5+0 



GXPF1 



6+, 1 

4+. 1 

5 + .0 



_6+, 1 
~4+, 1 



_6+, 1 
~4+, 1 



7+, 
0+ 1 



_6+, 1 
4+ 1 



6+ 1 
-4+ 1 



-2+ 1 

-1+0 

7+0 
~n+ i 



H u 



42 Ca 



2 Sc 42 Ti 



54 Co 54 Fe 



FIG. 2: (Color online) Energy spectra of two-particle states in the I/7/2 level predicted by the CD-Bonn (light blue), CD- 
Bonn+3terms (red) and GXPF1 (green) realistic interactions. Each is compared to the model Hamiltonian Hm (black) with 
X = 0.096, 0.124 and 0.092, respectively (Table \H\ . For comparison, the available experimental energy spectra of the A — 42 
Ca, Sc, Ti isobars (blue) and A — 54 Co and Fe isobars (magenta) are also shown. 



While correlation coefficients © prove useful in 
studies of nuclear properties shaped by the residual 
pure two-body interaction, the discrete energies of a 
quantum-mechanical system arc additionally influenced 
by the centroid, (fT2|) . and the overall interaction 

strength related to ([9]). The centroid of the Hm 

pairing+quadrupole model interaction, which coincides 

with the one for Sp(4) (W^ M = ) (Table [HJ, 

is very close to that for GXPF1 and both differ from the 
other two realistic interactions. However, this quantity 
is irrelevant for the energy spectra relative to the ground 
state within a given distribution. 

For the scalar distribution, the Hm model Hamiltonian 
generates energy spectra that are comparatively more 
spread out, especially with respect to CD-Bonn (see 
Figure [2] and a in Table HT|) . A rcnormalization of Hm 
will push the higher-lying states down and will establish 
an energy pattern very much like the ones observed in 
GXPF1, CD-Bonn+3terms and experiment (except for 
the 1 + state). Such a rcnormalization, however, is not 
done because the scalar distribution itself introduces av- 
eraging over isospin values. Instead we turn to the more 
detailed T = 1 spectral distribution. 

The T = 1 part of the pure two-body realistic interac- 
tions is reproduced quite well by the model Hamiltonian. 
This in turn yields a similar energy spectra (Figure [3]) 
as predicted by the model interaction Hm and both the 
CD-Bonn+3terms and GXPF1 realistic interactions. 
They also agree well with the available experimental data 
for the A — 42 and A = 54 isobars (the latter refer to 
systems of two holes). Here again, the agreement be- 
tween Hm and CD-Bonn is not as good as for the other 



two interactions, especially in reproducing the spread- 
ing of the states, which is smaller for CD-Bonn relative 
to the experimental energy spectra. Compared to the 
Sp(4) model interaction, the energy of the first 2 + state 
is brought slightly lower by the addition of Hq(2), to 
1.47 MeV (when \ is determined in comparison to CD- 
Bonn+3tcrms) and 1.22 MeV (to GXPF1). Such values 
are a bit closer to the experimental results, namely, 1.52- 
1.59 MeV for the A = 42 isobars and 1.41-1.44 MeV for 
the A = 54 isobars in the I/7/2 level. While the influence 
of the quadrupole-quadrupole interaction is significant in 
the 2 + T = 1 states, it does not affect the estimate for 
the Sp(4) parameters because they were determined with 
regard to the nuclear isobaric analog + states. 



<U 1 

a 
u 



CD-Bonn 



CD-Bonn+ 
3terms 



GXPF1 



lit. 



42 Ca 42 Sc 42 Ti 54 Co 54 Fe 



FIG. 3: (Color online) Energy spectra of T = 1 states pre- 
dicted by the CD-Bonn (light blue) , CD-Bonn+3terms (red) 
and GXPF1 (green) realistic interactions. Each is compared 
to the model Hamiltonian Hm (black) with \ = 0.071, 0.036 
and 0.055, respectively (Table IH) . For comparison, the ex- 
perimental T — 1 energy spectra of the A — 42 Ca, Sc, Ti 
isobars (blue) and A — 54 Co and Fe isobars (magenta) are 
also shown. 
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An interesting result is that in the detailed case of 
isospin-scalar distribution, the ( J = 0, T = 1) two-body 
matrix element is not affected at all by the Q ■ Q inter- 
action added to H sv ^ because Hq{2) has a zero con- 
tribution to this state. In this way, the addition of a 
quadrupole-quadrupole interaction does not alter the + 
T = 1 ground state of the A = 42 isobars, which are 
precisely the lowest isobaric analog + states where the 
Sp(4) model has been applied. Moreover, the average 
effect of this two-body collective interaction throughout 
the entire shell is zero [see (|19|)]. Even though this prop- 
erty of Hq{2) is imposed by construction, such a form 
(as already mentioned) is typically preferred so as to pre- 
serve the shell structure. In short, the domain of states 
where the Sp(4) model was applied is not influenced by 
the inclusion of the quadrupole degree of freedom to the 
pairing model. The Sp(4) model interaction itself actu- 
ally accounts for all the effects, small or large, due to 
the influence of the quadrupole-quadrupole interaction 
on these states. 



D. Correlations between Interactions for Nuclear 
Systems with More Than Two Nucleons 

An important feature of spectral distribution theory 
is that the correlation coefficient concept can be propa- 
gated beyond the defining two-nucleon system to deriva- 
tive systems with larger numbers of nucleons Q and in 
the isospin-scalar case, for higher values of isospin [l(|. 
The propagation formulae (|10|) determine how the aver- 
ages extracted from the two-nucleon matrix elements in 
the two-nucleon system get carried forward into many- 
nucleon systems. This propagation of information is 
model-independent. In this way one can track the sim- 
ilarity of pairing/rotational characteristics between dif- 
ferent interactions in many-nucleon systems (49j . 

In the scalar case the correlations between the inter- 
actions retain their values as given in Table U For the 
isospin-scalar distribution, the correlation coefficients be- 
tween the realistic interactions and the Hm model inter- 
action decrease for n > 2 and higher-T values compared 
to the n — 2 T = \ case and increase when compared to 
the n = 2T = estimates (Figure SJ). 




(a) (b) (o) 

FIG. 4: (Color online) Correlation coefficients, Ch' m h r > OI the 
pairing+quadrupole model interaction with (a) CD-Bonn, (b) 
CD-Bonn+3terms and (c) GXPF1 as a function of n and 
T for the I/7/2 orbit. The representation is symmetric with 
respect to the sign of n — 20,. 



For given n and T, the C values can again be found 
as the maximum correlation for an optimal value of the 
quadrupole-quadrupole strength Xi which is related to 
the angle between the geometrically represented Hr and 
its projection onto the Hm plane. When T = n/2 (Fig- 
ure 13} , which corresponds to the highest possible isospin 
states (including n = 2, T = 1) in 41 nuclei with va- 
lence nucleons occupying the I/7/2 orbit, the realistic 
interactions continue to be strongly correlated with the 

(T) 

Hm model Hamiltonian, namely, Ch h ^ s 0-95, 0.98 
and 0.96 for CD-Bonn, CD-Bonn+Sterms* and GXPF1, 
respectively. For these states, the other types of inter- 
actions are negligible and constitute only 3% of CD- 
Bonn+3terms, 8% of GXPF1 and 9% of CD-Bonn (Fig- 
ure QTb)). 

For all the cases with T 7^ n/2 throughout the en- 
tire shell, the correlation coefficients are found to retain 
almost the same value (Figure |4]), namely, Ch m Hr is 
around 0.63 - 0.70, 0.76 - 0.84 and 0.80 - 0.86 for CD- 
Bonn, CD-Bonn+3terms and GXPF1, respectively, with 
the corresponding optimal strength x of the quadrupole- 
quadrupole interaction in the intervals, 0.087 — 0.096, 
0.042 - 0.047, and 0.057 - 0.058. The smaller x, the 
weaker the Sp(4) symmetry breaking resulting from the 
additional quadrupole-quadrupole interaction. Again, 
the least strength is observed when Hm is compared to 
CD-Bonn+3terms. In the case of GXPF1, \ remains 
almost the same for all the states, with T 7^ n/2 as 
well as with T = n/2 (Table HH last row). In addi- 
tion, the relative contribution of other types of interac- 
tion, which are not accounted for by Hm f| 19[) . is some- 
what greater for CD-Bonn+3terms and CD-Bonn com- 
pared to GXPF1 because of their comparatively smaller 
Ch^h (Figure H}. In short, among the three realistic 
interactions, when T differs from n/2 the Sp(4) dynami- 
cal symmetry continues to be reflected the most in CD- 
Bonn+3terms, while the extended pairing+quadrupole 
model interaction correlates the best with GXPF1. 

An interesting result is that among the T ^ n/2 cases 
the highest H correlations for each realistic interac- 
tion are observed for the low-T mid-shell nuclear states, 
where the x strength is relatively smaller. This suggests 
that for these states the other kinds of interactions not 
present in Hm grow weaker and the realistic interac- 
tions are comparatively closer in behavior to the sym- 
plectic pairing H sv u\. The lowest correlation is observed 
for the mid-shell nuclei with isospin next to the highest, 
where x is the largest. For these cases, both quadrupole- 
quadrupole and other types of interactions accounted for 
by the realistic interactions increase in relative impor- 
tance. 

In summary, for all of the (n, T) distributions the 
model pairing+quadrupole interaction accounts on 
average for about 59%, 77% and 78% of the CD-Bonn, 
CD-Bonn+3terms and GXPF1 realistic interactions, 
respectively, and up to 91%, 97%, and 92% of those in 
the highest possible isospin states for a given n, where 
the Hm model interaction can be used to provide a 
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reliable description. 

Although, the aforementioned results refer to a single-j 
orbit (I/7/2), they represent a first step towards a gen- 
eralization to multi-j major shells (such as the fp shell 
or even a set of several major shells). Such an extension 
is certainly feasible and not space limited from the per- 
spective of spectral distribution theory. Even though the 
statistics for the I/7/2 level are not large compared to 
the fp shell, it is a natural choice. Specifically, the I/7/2 
orbital is comparatively far away from both the neigh- 
boring ds and upper fp shells which means there is a 
preponderance of a single-j coherence over configuration 
mixing in the low-lying nuclear states for such nuclei. For 
this reason it is an interesting example in its simplicity 
and provides a quite clear view of the pairing/rotational 
foundation of the nuclear interaction that is free of com- 
peting configuration mixing effects. In addition, it does 
represents a partitioning of the /p-space and as such it 
provides for more detailed spectral measures that may 
reflect important fine effects that are otherwise averaged 
out when the entire fp major shell is taken into account. 
It will be also important to augment this single-j shell 
study with a follow-on, complementary multi-shell anal- 
ysis of similar type. The multi-shell case will be the topic 
of a future publication. 

In addition to the above, we note for completeness that 
the estimates for the parameters of the model Sp(4) in- 
teraction are not really relevant to the primary objective 
of the present study. The reason is that the correlation 
coefficients in the detailed isospin-scalar case are inde- 
pendent of the interaction strength parameters, which 
therefore do not affect the correlation measures of the 
model interaction with realistic ones or the pairing and 
rotational characteristics of the latter. 



V. CONCLUSIONS 

With a view towards a broader study within multi-j 
shells, we compared three realistic interactions and two 
model pairing and quadrupole interactions for the I/7/2 
orbit by means of the theory of spectral distributions. 

In the more detailed case of isospin-scalar distribution, 
the CD-Bonn, CD-Bonn+3terms and GXPF1 realistic 
interactions were found to contain on average 59%, 77%, 
and 78%, respectively, of the pairing+quadrupole inter- 
action. Moreover, this percentage goes up to 91%, 97%, 
and 92%, respectively, for the highest possible isospin 
group of states for all the nuclei with valence protons 
and neutrons occupying the 1/7/2 shell. For these states, 



the strongest correlation was observed between the CD- 
Bonn+3terms and the pairing+quadrupole model inter- 
action, where other types of interaction accounted in the 
realistic interactions represent only 3% of it. They con- 
stitute 8% of the GXPF1 realistic interaction, and 9% 
of CD-Bonn. For these cases, the pairing+quadrupole 
model interaction has been shown to be a very good ap- 
proximation that provides a reasonable description of the 
energy spectra of the nuclei in the I/7/2 level. While both 
interactions, CD-Bonn+3terms and GXPF1, exhibit a 
well-developed pairing character compared to CD-Bonn, 
the latter appears to build up more (less) rotational col- 
lective features that are outside of the scope of the T = 1 
(T = 0) Sp(4) interaction. 

The major advantage of the sp(4) algebraic model, 
which focuses on the isovcctor pairing correlations and 
also includes a certain portion of the quadrupolc- 
quadrupole interaction, is that it provides an elegant so- 
lution for describing the pairing-governed isobaric analog 
+ states in light and medium mass nuclei. In addi- 
tion, while it correlates to a reasonably good extent with 
the realistic interactions, the description of the low- lying 
nuclear energy spectrum of higher- J states is improved 
with the inclusion of the quadrupole-quadrupolc interac- 
tion that being symmetric under SU(3) breaks the Sp(4) 
symmetry and removes degeneracies. Nevertheless, we 
found that for the isospin-scalar distribution the Sp(4) 
model interaction accounts for a large part of the CD- 
Bonn+3terms realistic interaction. It also includes be- 
tween 15% to 35% of the rotational collective interac- 
tion and typically accounts for a rather large portion of 
the overall correlation of the realistic interactions with 
the pairing+quadrupole interaction. Moreover, the addi- 
tional quadrupole degree of freedom incorporated in the 
symplectic model interaction does not affect the domain 
of states where the Sp(4) model was applied and hence 
introduces no errors in the estimates of the parameters 
of the symplectic interaction. These results confirm the 
conclusion that the Sp(4) interaction can provide for an 
approximate pattern of the nuclear energy spectra and, 
above all, can be accepted as a very reasonable approx- 
imation to describe the pairing-governed isobaric analog 
+ states in the nuclei in the I/7/2 orbit. 
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